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ABSTRACT 

We present galaxy circular velocity functions (GCVFs) for 34 low redshift {z < 0.15) 
clusters identified in the Sloan Digital Sky Survey (SDSS), for fifteen clusters drawn 
from dark matter simulations of hierarchical structure growth in a ACDM cosmology, 
and for ~29,000 SDSS field galaxies. We find that the observed and simulated cluster 
GCVFs take the form of a power law. The cumulative GCVFs of the simulated clusters 
are very similar across a wide range of cluster masses, provided individual subhalo 
circular velocities are scaled by the circular velocities of the parent cluster. Once 
all sources of error are accounted for, the intrinsic scatter in the cumulative, scaled 
observed cluster GCVF is consistent with the simulations. The slope of the observed 
cluster GCVF is '^-2.5, independent of cluster velocity dispersion. The average slope 
of the simulated GCVFs is somewhat steeper, although formally consistent given the 
errors. Using our highest resolution cluster, we find that the effects of baryons on 
galaxy rotation curves is to flatten the simulated cluster GCVF into better agreement 
with observations. Finally, we find that the field GCVF deviates significantly from a 
power law, being flatter than the cluster GCVF at circular velocities less than 200 
kms"-*^, and steeper at circular velocities greater than 200 kms~^. 

Key words: cosmology: theory — cosmology: observations — galaxies: halos — 
galaxies: clusters — methods: A^-body simulations — dark matter. 



1 INTRODUCTION 



In the widely accepted hierarchical clustering model of struc- 
ture formation, adiabatic fluctuations in the initial density 
field grow and merge into increasingly larger dark matter ha- 
los. Simple inflationary approaches model the initial power 
spectrum of fluctuations as Gaussian and scale-invariant, 
while cosmological parameters and the nature of the dark 
matter control its subsequent evolution. The resulting mass 
function of dark matter halos can be predicted using A^- 
body simulations JSheth fc Tormcn" 1999: 'Bode et alj|200lt 
Ijenkins et al.l l200lt iR^e^^t al. .2003j or analytic meth- 
ods based upon Press-Schecter theor y (|Press_^L.Schechten 
| l974l: 'Bond et al.' 'l99l'; 'Bowed Il99ll: iLacey fc Cold Il993l: 
ISheth . Mo & Tormcn 2001). The observed mass function is 
therefore an important constraint on fundamental proper- 
ties of the universe. In the following, we focus on the mass 
function on galaxy scales. 

Comparing predicted mass functions to observations re- 
quires mapping the observable properties of galaxies onto 
the masses of their dark matter halos. In one technique, 
semi-analytic models (SAMs) are used to predict the lu- 
minosities of galaxies within halos of different masses, by 
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parametrizing models of the complicated processes of gas 
cooling, star formation, stellar evolution, and feedback. 
These techniques efficiently explore a wide variety of mod- 
els. However, their utility is limited by the extent to which 
the simple recipes capture the relevant physics. 

An alternative to SAMs utilizes the empirical rela- 
tionship between velocities and luminosities in both spi- 
ral and elliptical galaxies to transform the observed galaxy 
luminosity function into a galaxy circular velocity func- 
tion (GCVF), which is closely related to the mass func- 
tion ( Cole fc Kaiser"l98^: 'Shima sakul[l993: IConzalez et alJ 
^000; Kochanck & White 2001). While this approach avoids 
a number of assumptions required to assign luminosities to 
dark matter halos, baryons cannot be completely ignored. 
The baryonic mass distribution affects a galaxy's circular 
velocity profile significantly, while also changing the equi- 
librium configuration of the host dark matter halo through 
adiabatic contraction iBlumenthal et al .Il98d) . Baryons also 
determine the probability that a galaxy will be included in 
a given observational sample. However, the circular velocity 
directly traces the total enclosed mass of a galaxy, which is 
dominated by dark matter at sufficiently large radii. There- 
fore, models of the CCVF are less sensitive to inaccuracies in 
our understanding of star formation, stellar evolution, and 
feedback than those of the luminosity function. The observed 
GCVF therefore better constrains the mass function. 
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Recent comparisons between observed and theoretical 
field GCVFs indicate that low-mass halos with circular 
velocities <, 2 00 kms ~^ are over-pre dicted by simulations 
jGonzalez etal,..200a : iKochanek fc W hite 200 L). It is now 
possible to extend the analysis of the GCVF to galaxies 
within clusters. The current Reneration of iV-body simu- 
lations Jchigna et alJliooi iTormen. piaferio fc Sveilll998l: 
i Klvpin et alJ Il999l: lOkamoto fc Habd I2OO0I: ISoringel et all 
12001!) has achieved the force and mass resolution necessary 
to capture the effects of tidal stripping, dynamical friction, 
and galaxy harassment , essential for predicting the GCVF 
in dense environments. iMoore et alj il999l) used the TuUy- 
Fisher relation to determine circular velocities for galaxies 
in the Virgo Cluster. They compared the resulting distri- 
bution to that in a standard CDM simulation of a Virgo- 
like cluster, and found good agreement down to ~50 kms~^ 
without accounting for the effects of baryons on galaxy rota- 
tion curves. Coupling dark matt er merging tree s with SAMs 
which include strong feedback, ISpringel et al.l ll200iri were 
able to reproduce a wide variety of cluster properties. Al- 
though they did not directly compute the cluster GCVF, 
the agreement between a variety of observed and modeled 
cluster relations (morphology-density, luminosity function, 
Tully-Fisher, Faber Jackson) suggests that they would have 
found agreement for the cluster GCVF as well. Although 
these works represent significant progress in testing simu- 
lations of cluster substructure against observations, larger 
samples of both simulations and observations are required 
to assess whether the results are typical, and over what clus- 
ter mass range. 



Toward this goal, we determine the GCVFs for 34 
nearby clusters from the Sloan Digital Sky Survey (SDSS). 
To compare the observed clusters to theoretical predictions, 
we assemble a sample of 15 high resolution cluster simula- 
tions, from which we measure the theoretical cluster GCVF. 
We also repeated the GCVF analysis for ~29,000 galaxies 
from the SDSS, to allow a comparison to the field. 



Our method for determining the GCVFs for the SDSS 
sample is described in In S|3|we describe the simulations. 
Our results are discussed in [21 First we discuss the shape 
of the cluster GCVF as related to that of the cluster lu- 
m inosity function (ijd .H. Second, we investigate the claim 
oflM oore et all (Il999l) that the velocity distribution of sub- 
structure on all mass scales is self-similar when scaled by the 
host halo velocity ( H4.2II . Third, we present power-law fits to 
the cluster GCVF, quantify its dependence on cluster veloc- 
ity dispersion, and make a comparison to dark matter clus- 
ter simulations ( >I4.3II . Fourth, we examine how the adiabatic 
contraction of dark matter halos due to baryons affects the 
comparison between observed and simulated GCVFs ( H4.4II . 
Finally, we compare the observed cluster GCVF to that mea- 
sured for the field ( H4.5II . A summary of our results can be 
found in 



Throughout we adopt Ho = 70 kms ^Mpc \ A = 0.7, 
and Q. = 0.3. 



2 METHODS: CONSTRUCTION OF THE 

OBSERVED CLUSTER GALAXY CIRCULAR 
VELOCITY FUNCTION 

To construct the observed cluster GCVF, we use existing 
cluster catalogs and data from the SDSS. For each cluster, 
we calculate the characteristic circular velocities for galaxies 
within the virial radius, using the Tully-Fisher and Funda- 
mental Plane relations. Note that this is done galaxy-by- 
galaxy, rather than by transforming the luminosity func- 
tion, as has been done previously ICo l e fc Kaiser. ■1981 
iGonzalez erai]l2000l: Ikochanek fc Whitj200li) In this sec- 
tion, we describe the SDSS, the cluster catalogs, our de- 
termination of global cluster properties, the application of 
the Tully-Fisher and Fundamental Plane relations, the com- 
pleteness velocities of the resulting GCVFS, and our back- 
ground subtraction procedure. 

2.1 The Sloan Digital Sky Survey 

The Sloan Digital Sky Survey (SDSS; lYork et a l."2000^ is 
a photometric and spectroscopic survey being carried out 
at the Apache Point Observatory usin g a 2.5-m telescop e 
equipped with modern CCD detectors llGunn et al.|[l998f) . 
This ambitious project will eventually cover the Northern 
Galactic Hemisphere, almost one quarter of the sky. The 
photometric survey provides ne arly simultaneous im aging 
in five passbands (it, g, r, i, z; iFukueita et"ailll99d) . The 
spectroscopic survey consists of medium resolution (A/ A A 
= 2000) spectra of galaxies, quasars, and stars selected sys- 
tematically from the photometric catalogs. 

The first public release of SDSS data occurred in June 
2001, and is known as the Early Data Release (EDR; 
IStoughton et alJl2002l) . The EDR consists of imaging data 
covering approximately 460 square degrees of sky, as well as 
approximately 55,000 spectra. 

The ph ot ometr ic and spectro scopic pipelines 
jLupton et alJ I2OOII: IStoughton et alJ ^002) measure 
numerous quantities that are available as part of the 
EDR. Below we describe only those which will be used in 
the present analysis. Note that because the photometric 
system has yet to be finalized, we denote the current, 
temporary magnitudes with asterisks, while references to 
the bandpasses remain simply italicized. 

The photometric pipeline computes two types of magni- 
tudes appropriate for use on galaxies: Petrosian and model. 
Petrosian magnitudes attempt to measure a constant frac- 
tion of the total light of a galaxy, independent of surface 
brightness, distance, or photometric model. The Petrosian 
radius is computed in the r band, and Petrosian magni- 
tudes in all bands are based upon the amount of fiux within 
two Petrosian radii. High signal-to-noise measurements of 
the Petrosian magnitude can be made for galaxies as faint 
as r* ~ 20. As Petrosian magnitudes become noisy, model 
magnitudes provide a more robust estimate of the total light 
of a galaxy. To compute model magnitudes, the photomet- 
ric pipeline fits a pure de Vaucouleurs profile and a pure 
exponential profile to the two dimensional r image of each 
galaxy. The best-fitting model is then applied, up to an am- 
plitude change, to images in the other bands. All fits account 
for the effects of local seeing. Parameters derived from these 
fits include r_deV, half the length of the major axis of the 
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ellipse that encloses half the light of the galaxy; and ab_exp 
and ab_deV, the best-fitting minor-to-major axis ratios for 
the two models. 

A bug in the SDSS photometric pipeline causes the 
model magnitudes of bright galaxies to be ~0.2 magni- 
tudes too bright. The inverse Fundamental Plane relation 
we use to estimate the circular velocities of early type galax- 
ies ( H2.4.21 was calibrated using model magnitudes afflicted 
by the same bug. We therefore continue in the use of EDR 
model magnitudes for early types in order to remain consis- 
tent with this calibration. The inverse TuUy-Fisher relation 
we apply to late type galaxies (12.4. H was calibrated using 
magnitudes extrapolated to infinity assuming exponential 
surface brightness profiles. Correctly-calculated model mag- 
nitudes would be the optimal choice for consistency with this 
calibration. In their absence, we use Petrosian magnitudes 
for late type galaxies. 

We adopt galacti c extinction corrections as function s 
of position based upon lSchlegel. Finkbeiner fc Davi3 l)l99c^ . 
and available from the SDSS database as the reddening 
parameter. All magnitudes in the remainder of this paper 
refer to reddening-corrected values. 

The spectroscopic pipeline determines both an emis- 
sion and cross-correlation redshift. Final redshifts (z) are 
assigned to each galaxy by choosing the method with the 
highest confidence interval. 

2.2 Cluster Sample 

Several cluster catalogs have been constructed from 
SDSS data, using a variety of identification algorithms 
jSahcall et ahllioOSl) . In order to compute accurate GCVFs 
down to the lowest possible circular velocities, we used the 
nearest clusters. Miller et al., 2002 (in prep) have identified 
62 clusters in the Equatorial Stripes of the SDSS usi ng the 
C4 algorithm ijNichol et al.ll200ll:lGomez et alJl2002h . This 
algorithm uses the observed phenomenon that galaxies in 
clusters have similar colors as well as locations. For each 
galaxy, the algorithm determines the number of neighbors 
in the seven-dimensional space spanning RA, DEC, redshift, 
and the four SDSS colors. This number of neighbors is com- 
pared to the distribution found for 100 random field galaxies. 
Based upon this comparison, the probability that the galaxy 
is a field galaxy is computed. Galaxies that are probably not 
field galaxies are identified using the False Discovery Rate 
thresholding technique (Miller ct al. 2001). 

In order to create as large a sample as possible, we sup- 
plement the C4 clusters with 51 nearby NED^ clusters in 
the same area. 

After removing 18 duplicates between the C4 and NED 
catalogs, 5 clusters with virial radii which extend beyond 
the SDSS data limits C ^Ol . 36 clusters (mostly NED) with 
fewer than 20 spectroscopically confirmed members ( >I2.3II . 
and 20 clusters with data defects or that were not rich 
enough to yield high signal-to-noise GCVFs, we are left with 
34 clusters, which are listed in Tabled 

1 The NASA/IPAC Extragalactic Database (NED) is operated 
by the Jet Propulsion Laboratory, California Institute of Tech- 
nology, under contract with the National Aeronautics and Space 
Administration. 



2.3 Determination of Global Cluster Properties 

As discussed in !^ the GCVF may depend on environment, 
and thus on cluster mass. While a reliable cluster mass esti- 
mate is difficult to obtain, the line-of-sight velocity disper- 
sion, (Tci, is relatively straightforward to measure, and can 
be used as a proxy for the mass. We calculate Ud using the 
following procedure. 

All galaxies with spectroscopic redshifts within a ra- 
dius of either 1.5 Mpc or the largest radius allowed by 
the extent of the data, whichever was smaller, were se- 
lected. Galaxies were ordered by recessional velocity, and 
gaps larger than 800 kms~^ were used to remove interlopers. 
A rob ust weighted gap estimator dBeers. Flvnn fc Gebhardl 
'1990?) was used to determine CTci from the remaining veloci- 
ties. This process was performed iteratively until the velocity 
dispersion stabilized. The resulting velocity dispersion was 
divide d by a factor of {1 -\- z) to correct it to the cluster 
frame (iDanese. de Zotti fc di Tullid[l98ol) . In order to en- 
sure accurate estimates of CTci, we rejected 36 clusters with 
fewer than 20 confirmed members. Errors were estimated us- 
ing bootstrap re-sampling, and are typically 10-20 per cent. 
The velocity dispersions and errors for each cluster can be 
found in Tabled columns 4 and 5, and the bottom panel of 
Figure |5| shows the distribution of CTci. 

In order to compare our circular velocity functions to 
each other and to models, we calculated them within a con- 
sistent aperture. The virial radius, i?vir, is a natural aperture 
to select. Girardi & Mczzctti ( 2003) present the following 
empirical relation between 7?vir and CTci: 

0.0017(7 /, -1 ,-1 N 
~ (^lj^z)'^/2 ^ ^wo Mpc). (1) 

Here Rvir is the radius within which the average density 
is equal to A ~ 200 times the critical density. Virial radii 
computed from Equation are listed in column 6 of Table 
Q We did not recalculate a^i using galaxies only within the 
virial radius for two reasons. First, for many of the clusters, 
nearly all spectroscopically confirmed cluster members were 
contained within -Rvir. Second, for the remainder of clus- 
ters, we wished to determine a^i using as many redshifts as 
possible. 



2.4 Determination of Galaxy Circular Velocities 

Each galaxy within Ryir and with r* < 21.5 (the star/galaxy 
separation limit; iLupton et al.l 120011) is given the desig- 
nation of Ear l y (E, SO, Sa) or Late (Sb, Sc, Irr) type. 
IStrateva et alJ (120011) find that u* — r* < 2.22 is optimal for 
separating these classes. For selecting early types, the com- 
pleteness is 68 per cent, while the reliability is 81 per cent. 
For selecting late types, the completeness is 77 per cent, and 
the reliability is 96 per cent. We then translated the surface 
photometry of each galaxy into a circular velocity using the 
conversions described in the following two subsections. 



2.4-1 Late Type Galaxies 

Disk galaxies exhibit a tight correlation between luminosity 
and circular velocity. The TuUy-Fisher relation (TFR) allows 
one to compute the former given the latter. We require the 
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Table 1. Galaxy Cluster Sample 
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RA and DEC are both in J2000.0; z is the cluster redshift; cr^i is the line-of-sight cluster velocity dispersion, as described in 
H2.3I AtTci is the bootstrap error on rj^i; Rvir is the virial radius of the cluster, as determined from Equation V^ompl is the 
completeness limit of the cluster GCVF, computed as described in il2.6l N2Q0 and /3 and their errors, AAf200 and A/3 are found by 
fitting Eauation ll6l to the cluster data; Af200 associated error, AA^jqq, is found by fitting Eg uation ll6l while fixing /3 to its 

average value of -2.5. 



inverse TFR, which converts from luminosity to circular ve- 
locity. Neither the TFR nor the inverse TFR h ave been cal- 
ibrate d in SDSS passbands. We therefore follow Isheth et all 
(2003) in adopting the inverse I-band relation presented by 
[GioyaneUi et al. Q992() for 360 spiral galaxies deemed to be 
members of 24 clusters (their "in" sample). Although it was 
calibrated using cluster data, we applied the following rela- 
tionship to late type galaxies in the determination of both 
the field and cluster GCVFs: 

log.o(2K)^2.5^-^- (^-^--^^|°g;o^-°°) , (2) 

where Mj is the A;-corrected, extinction-corrected I-band ab- 
solute magnitude, extrapolated to infinity assuming an ex- 
ponential surface brightness profile; and 14 is an estimate of 
the maximum rotational velocity of a galaxy. Although disk 
galaxy rotation curves exhibit a variety of shapes, the circu- 
lar velocity generally increases with radius until it reaches a 
broad peak at ~10 kpc, beyond which it remains relatively 
constant. Thus K should be similar to the circular velocity 
measured at the flat part of the rotation curve. 



In order to convert the observed SDSS magnitudes into 
those upon which the TFR is based, we first computed the 
r-band absolute total magnitude of each late type galaxy: 



Mr* = nir 



5logdL{z)+5-K{z), 



(3) 



where rrir* is the apparent r-band Petrosian magni- 
tude, d_t is the luminosity distance in parsecs, and K{z) 
is the fc-correction, computed using kcorrect v1_11 
^Blanton et al.ll2003h. We then use the conversion from 
iFukugita. Shimasaku fc Ichikawal ^1995^ for late type galax- 
ies: Mr" - Ml 0.9. 

Internal extinction corrections w ere estimated using the 
procedure of ITuUv fc Fouau3 dlS 



A 



-2.5 log 



/(l + e-") + (l-2/) 



1 - e" 



(4) 



where / = 0.1 is the fraction of stars homogeneously mixed 
with a dust layer having opacity r = 0.28; a = rseci; and i 
is the inclination, estimated as: 
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COS (l) = — 



(5) 



Here go = 0.20 is the intrinsic thickness and q — ab_exp 
is the r-band axis ratio. The extinction is considered to be 
constant for i > 80°. 

The internal extinction corrections were on the order of 
~0.5 magnitudes. 



2.4.2 Early Type Galaxies 

Early type galaxies occupy a thin plane in the space de- 
fined by stellar velocity dispersion, physical size, and sur- 
face brightness. Given the velocity dispersion and surface 
brightness, one can use the F undamental Plane rela tion (FP; 
iDiorgQvski fc Davis! (|l98i1) : IPressler etaP to esti- 

mate the size. In order to construct the GCVF from SDSS 
photometry, we require a relation that will yield the velocity 
dispersion (which is simply related to the circular velocity) 
given a size and surface brightnes s. 

In a series of three papers llBernardi et al.ll2003ll?ll?t) . 
Bernard! and collaborators constructed a sample of ~9000 
SDSS early type galaxies, from which they computed 
the joint distributions of tJios, the line-of-sight velocity 
dispersion measured through a circular aperture with a 
radius equal to one-eighth the half-light radius (ro = 
r.devV ab_dev) ; Jo, the average surface brightness within ro; 
and -Ro, the half-light radius in kpc. In order to determine 
the coefficients of the FP relation, iBernardi et alJ (l2003l) 
used the above joint distributions to minimize the sum of 
Ai = logj^Q _Ro — alogio cTios — blogio lo — c over all galaxies 
in their sample. In order to compute the coefficients for the 
inverse FP relation. 



logio cTios = a' log^o Ro + b' logjo 



Jo + c . 



(6) 



we used the same distributions to minimize the sum of 
AinvFP = (logio ""los ~ ^ logio -^0 - logio /o - c')^ Over all 
galaxies. We found a! = 0.597, h' = -0.447, and c' = 5.449. 

For each early type galaxy, wo estimated ctios using 
Equation |S| Following Bernardi et al. (200d), we set 



^0 = Mr. + 2.51og(27rr,1) - K{z) - 101og(l + z). 



(7) 



where [s.^ = -2.5 log lo and K(z) is as previously defined. 

In determining Ro, the physical characteristic size, we 
first measured the angular characteristic size for each galaxy, 
ro. Because the inverse FP that we used was calibrated 
against galaxies at a variety of redshifts, the angular sizes 
measured in each band were linearly interpolated in order to 
determine the angular size at a rest wavelength correspond- 
ing to the central wavelength of each filter. 

In order to plot early type galaxies on the same veloc- 
ity scale as late types, we converted aios to 14 (see H2.4.1II . 
Fitting non-parametric, spherical models to line profile and 
radial velocity dispersio n data for 21 elliptical galaxies, 
iKronawitter et all i2000f) found that the circular velocity 
profiles of elliptical galaxies are fiat t o the 10 per cent leve l 
outside a radius of 0.3 effective radii. iGerhard et alJ (120011) 
showed that there is a tight relation between the circular ve- 
locities in the flat regions of the rotation cu rves and the cen - 
tral velocity dispersions for these ellipticals. lFerrares3 (l2002|'l 
exam ined this relation, t aking central velocity dispersions 
from iDavies et al.l I^S^) for the ellipticals and correcting 



them to an aperture of one-eighth the effective radius to 
flnd a relation between logj^o ^"^^ log^o fios. For an isother- 
mal sphere, 14 = v^cios. Using the same data as iFerra resel 
(2002), we found that the isothermal assumption systemati- 
cally underpredicts 14, but that the data are consistent with 
the following relation: 



14 « 1.54aio 



(8) 



For each observed early type galaxy, we use Equation |H| to 
determine 14. 



2.^.3 Accounting for Scatter in the Derived Circular 
Velocities 

Values of 14 for individual galaxies were computed wit h the 
mean relations in t|2.4.2l and > l2.4.1llSheth et al.l ll2003l) note 
that if the scatter is significant, the use of mean relations 
can significantly distort the shape of the GCVF. Therefore, 
in constructing our GCVFs, we accounted for the scatter 
by representing each galaxy's 14 as a probability distribu- 
tion, rather than a single value. Because we calculated 14 
using different relations for early and late type galaxies, this 
distribution depends on galaxy type. 

The circular velocities for late type galaxies were com- 
puted directly from the inverse TFR ( t|2.4.1ll . We therefore 
adopted the model fo r the scatter around 14 provided by 
iGiovanelli et al.1 (ll997^ . 

2 ^_ (-0.28(logio(2K)- 2.5) +0.26)" 
< Ai„.TF >- . (9) 

The circular velocities for early type galaxies were com- 
puted using Equations El and IHl There were two sources of 
scatter in determining 14: the scatter around aios and the 
uncertainty in the relationship between uios and 14. We esti- 
mated the former by excha nging V and R i n the expression 
for < Al > in Equation 4 of lBernardi et all |2003) and using 
the joint distributions of crioa, Ro, and lo presented in their 
Table 1. Doi ng so yielded a sca tter of < Af„^pp 0.061 
in logj^Q (Tios. ISheth et al.l ll2003l) found that the scatter in 
logiQ (Tios depends slightly on log^g trios. This was a small ef- 
fect compared to our other sources of error; we therefore 
treated the scatter in the Fundamental Plane as constant. 
The conversion from crios to 14 (Equation |HJ was calibrated 
on only ~20 elliptical galaxies. Therefore a robust model of 
the scatter due to this conversion was not possible and we 
neglected it. 

Because the Fundamental Plane and TuUy-Fisher re- 
lations were calibrated using samples created under more 
stringent selection criteria than we used, the above relations 
underestimate the scatter that exists in our data. 



2.5 Background Subtraction 

The largest uncertainty in our determinations of cluster 
GCVFs stemmed from the difficulty in determining clus- 
ter membership. Because a spectrum is not available for 
every galaxy in the field of each cluster, our background 
subtraction is statistical. We calculated the circular veloc- 
ity for every galaxy within the virial radius of each clus- 
ter under the assumption that it lies at the cluster redshift. 
Clearly, some galaxies have redshifts greater or less than that 
of the cluster, and our calculations were incorrect for these 
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galaxies. This effect was removed during our background 
subtraction procedure. Again using the cluster redshift, we 
calculated the circular velocity for every galaxy within an 
annulus bounded by inner and outer radii of 3 and 5 Mpc, 
respectively. These galaxies should have a redshift distri- 
bution similar to that of the cluster interlopers. Thus, by 
subtracting the appropriately normalized GCVF of galaxies 
in the annulus, we corrected for the interlopers. 

The error in the number of galaxies in each circular 
velocity bin are calculated assuming Poisson noise in both 
the total counts within the virial radius, Ctot , and the counts 
within the background annulus, Cbackground: 

Cat = Ctot + / Cbackground, (10) 

where / is the normalization factor applied to the back- 
ground annulus counts to account for the difference in 
its area and that of the virialized part of the cluster 
(0.02 < / < 0.3). 

The contribution to the total error in the GCVFs from 
the Poisson error in the background counts decreases with 
increasing Vc- Above ~200 kms~^, the contribution is less 
than 40 per cent. 

2.6 Completeness 

Because our procedure distinguishes between morphological 
types, we could assign different values of T4 to galaxies of the 
same luminosity. Our requirement that r* < 21.5 therefore 
leads to a cut in Vc that depends on both redshift and galaxy 
type. Figure illustrates the effect our magnitude limit has 
on Vc for a representative cluster. The filled points represent 
early type galaxies, while the hollow points represent late 
type galaxies. For this cluster, 90 per cent of galaxies with 
21 < r* < 21.5 have Vc > 83kms~^. Below this completeness 
velocity, V^ompi, significant numbers of galaxies are missing 
from our sample. All GCVFs in what follows are truncated 
below Kompi, which was calculated independently for each 
cluster. The adopted values are listed in Table Q column 7. 



3 METHODS: SIMULATIONS 

In order to compare our observational results against simula- 
tions, we compiled a set of 15 simulations of the dark matter 
components of clusters. Each simulated cluster was evolved 
in the currently favored ACDM cosmology with A=0.7, Ho 
— 70 kms~^ Mpc~^, flo = 0.3, and erg = 1. Ten are the most 
massive clusters from the zero redshift output of the 50^ h^^g 
Mpc'^ simulation described in lReed et al.l ^l2003^ . This simu- 
lation has a particle mass of 1.9x IO^Mq and a force resolu- 
tion of 7.14 kpc. An additional four clusters were simulated 
using the volume renormalization technique jKatz fc White! 
[l993) and have (0.5-0.9) X W particles within the virial ra- 
dius and force softenings of <, 1 per cent of the virial radius. 
Further detail s on t hese four simulations can be found in 
iBoreani et all tOQ'j) . 

The 14 cluster simulations were of sufficient resolution 
to robustly determine peak circular velocities in subhalos 
above the completeness velocities of our observed GCVFs 
(see >12.til . They were of insufficient resolution, however, 
to accurately trace subhalo internal density profiles. These 
were necessary to determine the effects of baryons on galaxy 
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Figure 1. Illustration of the completeness limit for one cluster 
in our sample. Each point represents a galaxy within R^ir for a 
representative cluster in our SDSS sample. No background sub- 
traction has been carried out. Filled circles are early type galaxies, 
and hollow circles denote late type galaxies. Ninety per cent of all 
galaxies with 21 > r* > 21.5 lie below the completeness velocity 
of Vcompl = 83 kms^^ for this cluster. 



Table 2. Simulation Sample 





0"cl 

kms~^ 


M200 
IO^^^Mq 


R200 
/i7oMpc 


10* Mq 


10'^ 


e 

kpc 


BoxOl 


525 


29.4 


1.73 


1.9 


1.5 


7.14 


Box02 


458 


26.6 


1.67 


1.9 


1.4 


7.14 


Box03 


383 


21.7 


1.56 


1.9 


1.1 


7.14 


Box04 


456 


18.4 


1.48 


1.9 


0.97 


7.14 


Box05 


451 


17.9 


1.46 


1.9 


0.94 


7.14 


Box06 


342 


15.4 


1.39 


1.9 


0.81 


7.14 


Box07 


388 


13.0 


1.32 


1.9 


0.68 


7.14 


Box08 


332 


11.6 


1.27 


1.9 


0.61 


7.14 


BoxOS 


369 


10.3 


1.22 


1.9 


0.54 


7.14 


BoxlO 


366 


8.9 


1.26 


1.9 


0.47 


7.14 


Hickson 


234 


2.5 


0.76 


4.41 


0.0486 


2.5 


Fornax 


400 


5.9 


0.65 


4.41 


0.118 


2.5 


Virgo 


565 


30.4 


1.75 


14.9 


0.179 


1.25 


Coma 


1007 


133.6 


2.86 


119.2 


0.090 


3.75 


HR Virgo 


600 


30 


1.75 


3 


1.00 


1.25 



Further details on the simulated clu sters named BoxOl - BoxlO 
can be found in iReed et al.l j20o3). The simul ations labelled 
Hickso n, Fornax, Virgo, and Coma are described in lBorgani et all 
f'2Q02|) . The last is used in i|4.4l Columns: ad is the velocity 
dispersion of the simulated cluster, -R200 is the radius within 
which the average density is 200 times the critical density, M200 
is the mass within -R20O1 i^DM is the mass of each dark matter 
particle, is the number of dark matter particles within 

i?200; and e is the force softening. 



The Cluster Galaxy Circular Velocity Function 7 



Simulated Clusters 



SDSS Clusters 



200 400 600 800 1000 1301 

c^oi (km/s) 

Figure 2. Top Panel: The distribution of cluster velocity disper- 
sions for the fifteen simulations described in J21 Bottom Panel 
The distribution of cluster velocity dispersions, calculated by the 
methods described in >I2.3I for our sample of 34 SDSS clusters. 



rotation curves, as we describe in H4.4I For this purpose, we 
used a very high resolution A'-body simulation of a galaxy 
cluster with CTcI = 600 kms^^ and Mvir = 3 x IO^^Mq, rep- 
resented by approximately one million particles within _Rvir 
= 1.75 Mpc. It was evolved in a 100'^ Mpc^ volume where 
the force spline softening was 1.25 kpc and the large scale 
structure was sampled at a lower resolution. 

Subhalos were id entified using SKID JStadeJ I2001I: 
iGovernato et al]ll997^ . a halo finder based on local density 
maxima and particulary suited to finding subhalos within 
larger structures. Peak circular velocities were also output 
by SKID. We used Tipsy^ to extract full density profiles 
from our highest resolution simulation. 



4 RESULTS 

4.1 The Shape of the GCVF 

A Press-Schechter mass function combined with M oc V'^ 
and a value of n = 2 for the fiuctuation-spectrum index give 



oc K 



pk 



dN (VpQ 

iMoore et all lll99 



(11) 



and lGhigna et al.l J2000h found the 
GCVF in high resolution simulations of a Vi rgo-like clus- 
ter t o be consistent with the above power law. iMoore et alJ 
il999,l also found the observed GCVF of Virgo to be consis- 
tent with it. In Figure 13 we plot the GCVFs for our sample 
of 34 clusters in three cluster velocity dispersion bins. The 
bins were chosen such that each contains approximately the 
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= -2.784±0.060 
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Figure 3. The observed CGVF, binned by cluster velocity disper- 
sion. Top panels: Individual cluster GCVFs are shown with thin 
lines. The dark line in each panel has a slope of -4, and is provided 
for reference. Cluster velocity dispersion ranges are shown in the 
upper right of each panel, in kms~^. Bottom panels: Each panel 
shows the composite functions constructed from clusters having 
cF^l in the range indicated in the panel directly above. Error bars 
take into account Poisson error in counting galaxies both within 
iJvir and in the background annulus. A power law fEauation llfil 
was fitted to these data, and is shown as a thin solid line. The 
slopes of the composite functions are shown in the upper right 
hand corner of each panel. 



same number of clusters. The bold lines have a slope of -4 
and are for reference only. From this figure, we can see that 
the observed cluster GCVF is a power law. The power law 
form of the cluster GCVF, while predicted analytically and 
numerically, is nonetheless somewhat surprising. The clus- 
ter luminosity function is usually described by a Schechter 
function. Because the luminosities of both early and late 
type galaxies are related to their circular velocities by power 
laws {L oc V^), it is not obvious that the galaxy population 
as a whole should obey both a Schecter luminosity function 
and a power law GCVF. 

Key to the explanation is the fact that both early and 
late types contribute to the GCVF, and that they have both 
different relations between Mr* and Vc, and make up dif- 
ferent fractions of the population with magnitude. Because 
the slope of the 14 — Mr relation for early types is steeper 
than for late types, late type galaxies were mapped onto a 
smaller range of velocities than early types. Thus, late type 
galaxies that populate the knee in the luminosity function 
were mapped to lower velocities than early type galaxies 
with the same luminosity. The proportion of early and late 
type galaxies that make up the knee in the luminosity func- 
tion affects the placement and prominence of the knee in 
the GCVF. In general, it was smeared out in the GCVF. 
This is the main reason that our observed cluster GCVFs 



8 



V. Desai et al. 




a 2.5 3 

log.o (Ve/(km s-i)) 

Figure 4. The composite CGVF, constructed from all 34 SDSS 
clusters in our sample. Top panel: The error bars were computed 
as m FigureEl The best-fitting power law is shown as a thin solid 
line, and its slope is shown in the upper right hand corner. Bottom 
panel: The residuals of the composite cluster GCVF from the best 
fit power law. The power law is a very good fit to the composite 
GCVF, except for the lowest Vc bin. 



look like power laws. If a single luminosity-velocity relation 
is used to compute a GCVF from a population of galaxies 
drawn from a Schechter funct ion, the resulting GC VF will 
also be a Schechter function JCole fc Kaiseilll989h . If two 
diflerent luminosity-velocity relations are used in the com- 
putation, it is still possible to arrive at a Schechter GCVF, 
as fou nd by Gonzalez ct al. (20o3) and [Kochanck & White 
i200l[) . A Schechter GCVF will not necessarily result, how- 
ever, depending upon the relative abundances of early and 
late type galaxies as a function of luminosity. 

Another reason for the power law GCVF is that our 
completeness limits did not allow us to observe the shape of 
the GCVF at velocities dominated solely by galaxies from 
the flat part of the luminosity function, where we would 
expect the GCVF to flatten out. Thus, very few bins define 
the knee. 

Finally, the signal-to-noise in each cluster GCVF is rela- 
tively low due to the small number of galaxies used. It would 
be difficult to trace the high velocity exponential drop-off if 
the GCVF did follow a Schechter function. 

In order to increase the signal to noise in in our de- 
termination of the cluster GCVF, we built composi te clus- 
ter G CVFs using a technique analogous to that of ICoUesl 
l|l989fl for building composite cluster luminosity functions. 
The number of galaxies in the jth Vc bin of the composite 
cluster GCVF is given by: 



_ iV,o V- iV»i 



where Nij is the number of galaxies in the jth bin of the 
ith cluster GCVF, TV^o is the normalization used for the 
ith cluster GCVF (taken as the number of galaxies brighter 
than 100 km s"^), rrij is the number of clusters contributing 
to the j'th bin, and A'^co is the sum of all the normalizations: 



The error on Nrs is: 



Ae, 



El 



/AA,, 



(13) 



(14) 



where /S.Nij is the error on the jth bin of the ith cluster. 

The composite cluster GCVFs in the three o"ci bins are 
shown in the bottom panel of Figure |3 The composite clus- 
ter GCVF constructed using all 34 clusters is shown in Fig- 
ure^] An increase in signal to noise still yielded a power law 
form. 



4.2 Is cluster substructure self-similar? 



iMoore et all Jl999t) found that the cumulative GCVF for a 
Virgo-like simulation is identical to that of subhalos within 
a Milky- Way-like simulation, provided the velocities are 
scaled by the velocity of the parent halo {Vd ~ \/2crci). 
In the top panel of Figure |5] we have plotted the cumula- 
tive, scaled GCVF for all fifteen simulated clusters from Ta- 
ble |5] Our simulated clusters have very simila r cumulative, 
scaled GCVFs, in agreement with the results of lMoore et alJ 

In the bottom panel of Figure|^ we plot the same quan- 
tity for our sample of SDSS clusters. The slopes and offsets 
of the simulated and observed clusters cannot be compared 
directly, since the effects of baryons have not been included 
in the simulations. However, assuming that baryonic physics 
changes the shapes of all the GCVFs in a similar way, we can 
compare the scatter in the observed and simulated GCVFs. 
The observations display a significantly larger scatter than 
the simulations. Uncertainties in the measurement of Ki and 
Poisson errors both contribute to the observed scatter. Tak- 
ing the plot of logj^Q[A(> (K/Ki))] versus Vc/Vd to be a 
power law, the variance in log]^Q[A(> (K/14i))] due to mea- 
surement errors in CTcI is 



1 A 



Mogio[iV(>(Vc/Ki))l 



In 10 



(15) 



(12) 



where Vd = \/2aci- Using the values of Ud found in 
and computing Av„, from the errors in iTci, we found that 
the scatter due to measurement errors in ad is 0.009, while 
the scatter due to Poisson noise is 0.012. The total observed 
scatter is 0.064. Thus, the instrinsic scatter is ~0.043. This is 
very similar to the intrinsic scatter we measured in the simu- 
lations. Thus our observed GCVFs show the same degree of 
variations found in ACDM dark matter cluster simulations. 



4.3 Trends in the Cluster GCVF with Cluster 
Velocity Dispersion 

Figure |3| shows that the cluster GCVF depends upon the 
cluster velocity dispersion. To quantify the observed cluster 
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Figure 5. The simulated and observed SDSS CGVFs, with 
galaxy circular velocities scaled by cluster circular velocites. The 
error bars have length 2A, where A is the RMS scatter at 
l°gio(^c/Vci ~ —0.7) due to various sources, as labelled. Pois- 
son refers to Poisson noise in counting galaxies within i?vir for 
the simulated clusters and to Poisson noise in counting galaxies 
within i?vir and the background annulus for the SDSS clusters. 
Measured refers to the total measured scatter. Al^l refers to the 
scatter in log^Q[Af(> (Vc/Vd))] due to errors in measuring Vd, 
and Intrinsic refers to the portion of the total measured scatter 
unaccounted for by Poisson scatter in the case of the simulations 
or by both Poisson scatter and that induced by AVci for the SDSS 
clusters. 



GCVF and its dependence on a^i, we fitted each with a 
function of the following form; 



log 



diV(K) 



loe 



iV2oo + /31ogio ( aoo^ms-O ' ^^^^ ^^oo = ( 



where (5 represents the slope of the power law and 7V200 
is the normalization at 14=200 kms~^. The results of the 
fit are reported in Table Q columns 8-11, and are plotted 
against logj^QO-c! in Figure |S| (filled points). The top panel 
of Figure |S| shows that the slope, /3, while showing a fair 
amount of scatter, does not vary systematically with 0^1- 
It has a mean value of /? = —2.5 ± 0.8. The slope of the 
composite cluster GCVF shown in Figure|l|is —2.61 ± 0.04, 
consistent with the mean slope. The normalizations of the 
individual cluster GCVFs are plotted in the middle panel of 
Figure|S|and show a clear systematic trend in the sense that 
higher mass clusters host more galaxies. The best fit line to 
these data are 



logjo iVaoo = (2.26 ± 0.2) log^Q - (6.84 ± 0.6). 



(17) 



Finally, we recalculated our normalizations with the slope 
fixed to its mean value of —2.5. These values of A^200 can be 
found in TableQ columns 12 and 13, and are plotted against 
(Tci in the third panel of Figure |S| The best fit line to these 
data are 




2.6 2.8 

l°gio('^ci/(km s-')) 

Figure 6. Each observed and simulated cluster GCVF 
was fit to the equation log-^Q dN / dVc = logj^g Af200 + 
/31ogio (14/200 kms^i). Top Panel: The slope of the cluster 
GCVF as a function of cluster velocity dispersion. Solid points are 
for the observed cluster GCVFs. The mean value is indicated by 
the solid line. Open points are for the simulated cluster GCVFs. 
The mean value is indicated by a dashed line. Middle Panel: The 
normalization of the observed cluster GCVF as a function of clus- 
ter velocity dispersion. High a^i clusters tend to contain more 
galaxies than low a^i. The line that best fits this trend is given 
by Eauation ll7l and is overplotted. Bottom Panel: The normaliza- 
tions of the fits for the observed cluster GCVFs were recomputed 
while holding the slope fixed to the mean value shown in the top 
panel. The new normalizations are plotted against a^i in the bot- 
tom panel. The line that best fits these data is given by Equation 
1181 and is overplotted. 



2.42 ±0.2) o-ci - (7.37 ±0.5). 



(18) 



This slope is somewhat smaller than theoretically expected, 
but not at a statistically significant level. 

We compared the slopes of the observed GCVFs to 
those of the simulated GCVFs. The latter are plotted as 
open circles in the top panel of Figure |S] The dashed line is 
drawn at their mean value: Paim = —3.4 ± 0.8. The average 
value of the simulations is somewhat steeper than that of the 
observations, although formally they agree within the uncer- 
tainties. It should be noted that the slopes are not directly 
comparable because the simulations do not account for the 
effects of baryons. The effects of baryons on the slope of the 
GCVF are discussed in the next subsection. 



4.4 The Effects of Baryons on the Cluster GCVF 

We found in i|4.3l that subhalos in simulated clusters have 
steeper GCVFs than measured for galaxies in clusters. We 
also noted that a more meaningful comparison between ob- 
servations and simulations would require a determination of 
the circular velocity profile taking into account the effects 
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of both dark matter and baryons. Baryons affect the cir- 
cular velocity profile of a galaxy in two ways. First, they 
significantly contribute to the enclosed mass at small radii. 
Second, they adiabatically pull the dark matter into a more 
concentrated equilibrium configuration than it would have 
in their absence. The total circular velocity profile is there- 
fore given by 

K'W = V;%(r) + K'DM(r), (19) 

where Vc,b(7') is the circular velocity profile of the collapsed 
baryonic mass and Vc,uM(r) is the circular velocity profile 
of the adiabatically contracted dark matter halo. 

The degree of these effects depends on the mass distri- 
bution of the baryons. The challenge lies in populating the 
subhalos with galaxies in a way that mimics, as closely as 
possible, observed empirical relations. A full description of 
the baryonic mass distribution requires knowledge not only 
the density profile shape (galaxy type), but also a normaliza- 
tion (total galaxy mass) and a scalelength (galaxy size) . Our 
method for populating the subhalos identified in our highest 
resolution simulation is described in Appendix^ Given the 
density profiles of the baryonic component of each subhalo, 
their circular velocity profiles were computed. Hernquist 
profiles are spherical, so the rotation curves of early type 
galaxies were straightforward to calculate: Kph = GM(r)/r. 
Determining the rotation curve of exponential profiles was 
also staightforward, but required taking account of the non- 
spherical geometry. 

The determination of Vcomir) required modelling how 
the dark matter responds to the collapse of the baryons. 
We assumed that the baryons contract slowly and that 
the halo remains spherical. Under these conditions, a dark 
matter particle conserves angular momentum as it moves 
closer to the center of the potential ('Blumenth al et al]ll986l : 
[Dalcanton et al. 1997; Mo. Mao fc White 199al : 

Mf(rf)rf = Mi(n)ri, (20) 

where ri is the initial mean radius of the particle, ri is the 
final mean radius of the particle, and 

Mf(rf) = Afb(rf) + M(n)(l - /cod). (21) 

Equations I2UI and 1211 both require knowledge of M(ri), the 
initial mass profile of the dark matter halo, before adiabatic 
contraction. While Vpk was robustly measured from all of 
the simulations listed in Table |21 a reliable determination 
of V{ri) or M (ri) could only be made from the simulation 
described by the last entry in Table |5| which has both very 
high resolution and a small softening scale. Using the mea- 
sured M(ri), the values of /cooi described in i|A0.2l and the 
Mb(r) computed from the density profiles described in Ap- 
pendix^ we calculated M(rf) from Equations I2UI and 1211 
From M{jf) it was trivial to calculate 14,dm(''). 

Once the total circular velocity profile Vdj-) due to dark 
matter and baryons was determined, we measured a char- 
acteristic circular velocity, K, that is analogous to that de- 
duced for our observed early and late type galaxies. We es- 
timated Vc as the circular velocity at three scale lengths. 

In Figure |7| we compare the simulated GCVF to that 
of SDSS clusters that have similar velocity dispersions. The 
dashed line represents the simulated cluster GCVF in the 
limit of no baryons and assuming that Vc = Vpk. The simu- 
lated cluster GCVF in this case is too steep, underpredict- 



ing the number of galaxies at high Vc. The bold solid lines 
bracket simulated GCVFs which have been corrected for the 
effects of baryonic infall. We varied the total mass that in- 
falls to form ellipticals (0.025 < /cool < 0.1) and the peak 
circular velocity above which subhalos host early type rather 
than late type galaxies (125kms~^ < Vbroak < 150kms~^). 
See Appendix^ When the effects of baryons were included, 
the number of high Vc galaxies increased, while the num- 
ber of low Vc galaxies decreased, bringing the simulations 
into rough agreement with the observations. This effect held 
whether Vc was measured at 2, 3, or 4 scale lengths, but was 
strongest if V was measured at 2 scale lengths. 

While ACDM models can reproduce the cluster GCVF, 
they overprodu ce the number of field galaxies with 
V< 1 20 kms~^ llGonzalez et al.l l2000t iKochanek fc Whitd 
l2001f) . Any solution to the disagreement in the field 
must preserve the agreement found in clusters. Modi- 
fications to CDM would affect the substructure within 
clusters as well as in the field, and tend to contra- 



diet other observations i 
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Supernova feedback alone fails to reproduce the luminosity 
function of the Local Group llSomervillel2002l) . but a combi- 
nation of feedback and squelching can. However, squelching 
cannot play a role at the circular velocities we have inves- 
tigated here. It affects only galaxies with circular veloci- 
ties <, 50 kms~^. The importance of supernova feedback 
at higher circular velociti es is uncertain. Recent work by 
iMac Low fc Ferraral (1199911 suggests that it can only regulate 
star formation, not remove mass at these scales. Regardless, 
it should be equally important in clusters and in the field. 
Finally, any incompleteness that exists in our cluster GCVF 
should also exist in our field GCVF. 

4.5 The Cluster Versus Field GCVF 

The circular velocity function of dark matter halos within 
clusters is expected to be a function of environment. Al- 
though we found only small differences among the pop- 
ulations of massive clusters, there may still be substan- 
tial differences between clsuters and the field. First, hier- 
archical clustering theory predicts that massive halos form 
from rarer density peaks and should therefore be more 
cluste red than low mass hal os (Kaiser 1984; Davi s et alJ 
1198^ IWhite et allllQSl Il98a: .Governato et al...l99^ . Sec- 
ond, halos within clusters are subject to the mass-altering 
processes of tidal stripp ing, dyn amical friction, and galaxy 
harrassment (Mooro et al][l99y). Third, low mass halos in 
dense regions formed, on average, earlier than those in the 
field. Reionization can therfore preferentially suppress the 
late-time formation of galaxies within low mass halos in 
the field fOuinn. Katz fc Efstathiou"l99el: ItuIIv et aljliooi : 
Benson et al. 2002; Somervillc 2002), although at velocites 
belo w what we can probe in this work. 

ICole fc Kaiseil il989r) determined the circular velcoity 
function of field galaxies by starting with a Schechter lumi- 
nosity function and assuming L oc V^~^, not taking into 
account that a late type galaxy has a lower circular ve- 
locity than an early type with the same absolute magni- 
tude, or that their relati ve abundances depend upon lumi- 
nosity, jshimasak^ (ll993;) computed the field GCVF galaxy- 
by-galaxy using both HI 21 cm line widths (~400) and the 
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Figure 7. Comparison of the GCVF measured from our highest 
resolution galaxy cluster simulation (HR Virgo) to that measured 
in SDSS clusters. The dashed line is the simulated GCVF con- 
structed using the Vpi; values output by SKID. Thick solid lines 
represent the range of corrected GCVFs given uncertainties in 
how gas infall affects galaxy rotation curves (see Appendix IXI for 
details). Thin solid lines are the GCVFs of SDSS clusters with 
measured velocity dispersions within 25 kms~^ of that of the 
simulated cluster. 



Tully-Fisher relation (-^300) for spiral galaxies and stellar 
velocity dispersions (~100) and the Faber Jackson relation 
(~90) for early type galaxies. They found that early types 
dominate at the high veloc ity end. 

.Gonzalez et al.l (l200d) determined the field GCVF by 
performing variable transformations on the type-specific B- 
band luminosity functions determined from the SSRS2 sur- 
vey. They found that, given the available Tully-Fisher and 
Faber Jackson Relations, all field galaxies could be treated 
as late types without changing the results. Doing so, they 
found that the number density of galaxies per unit velocity 
can be described by a Schecter function. 



*(Vc)dK = ** ( I exp 



ii 



where 4-. = (3.2 ± 0.6) x 10"^ Mpc~^ = -1.3 ± 0.13, 
n = 2.5, K,* = 247 ± 7. iKochanek fc White! J200l[) found 
similar results, although they took into account the di ffer- 
ent tr ansformations for early and late types. Sheth et alj 
i2003h analyzed field early type galaxies within the SDSS, 
and found that accounting for the scatter around the mean 
Tully-Fisher and Faber Jackson or Fundamental Plane rela- 
tions is extremely important for reproducing the high veloc- 
ity end of the GCVF. 

In order to allow a direct comparison to our cluster 
GCVFs, we constructed a field GCVF using the same rela- 
tions that we used in the previous cluster analysis. We se- 
lect galaxies from the SDSS EDR spectroscopic sample with 



14.5 > r* > 17.77. Each of the resuhing ~29,000 galaxies 
was assigned a weight equal to the inverse of the volume 
over which the galaxy could be observed, given the appar- 
ent magnitude limits. Circular velocities were computed for 
each galaxy as described in >I2.4.1I and H2.4.2I The scatter 
was dealt with as described in H2.4.3I The value of the field 
GCVF in a given circular velocity bin is equal to the sum 
of the weights of the galaxies in that bin. The results are 
plotted as points with error bars in Figure El along with the 
functi onal forms for the field GCVF given bv lConzalez et aD 
hood) (Equation mi solid bold line). Figure [H illustrates 
that the cluster GCVF shows much less curvature than that 
of the field, being steeper for Vc ^ 200km s"'^ and more shal- 
low for K > 200 km s-^ 

Following the discussion in Section r4.3l the difference in 
the shapes of the field and cluster GCVFs can be explained 
by their different relative abundances of early and late type 
galaxies as a function of luminosity. It must be noted that 
our division of galaxies into early and late types depends 
upon the u* — r* color, which is sensitive to star forma- 
tion. There is some evidence that the cluster environment 
suprcsscs star form ation (Couch et al. 1987; Bargcr ct al\ 
\l99dL IPoggianti et^ alJll999l) . A field galaxy that falls into 
a cluster and experiences truncated star formation would 
redden, but its r*-band luminosity would not change signif- 
icantly, because r* is sensitive to the total stellar mass. If 
the change in u* — r* is large enough to change late types 
into early types, the cluster GCVF would have more high- 
velocity galaxies and fewer lower- velocity galaxies relative to 
the field, as is observed. While the trend is correct, detailed 
models are required to determine to what extent this effect 
can explain the observed differences between the field and 
cluster GCVFs. 



5 SUMMARY 

We constructed the cluster and field galaxy circular veloc- 
ity functions using data from the Sloan Digital Sky Sur- 
vey Early Data Release. Circular velocities were determined 
using the inverse Tully-Fisher and Fundamental Plane rela- 
tions on a galaxy-by-galaxy basis, rather than by performing 
a transformation of variables on a functional fit to the lumi- 
nosity function. For comparison to the observations, we also 
analyzed the GCVFs determined from 15 clusters simulated 
in the ACDM concordance cosmology. Our results can be 
summarized as follows: 

1. We found that, although the observed cluster lumi- 
nosity function is a Schechter function, the observed cluster 
GCVF can be described as a power law. 

2 . We tested the self-similarity claimed bv lMoore et alJ 
il999() for our scaled simulated and observed cluster GCVFs. 
The simulated GCVFs are indeed very similar to one an- 
other, and show only a small scatter. The scatter in the ob- 
served cluster GCVFs is much larger, but can be explained 
by measurement error. Thus, the intrinsic scatter in the ob- 
served cluster GCVF is similar to that predicted. 

3. The normalization of the observed cluster GCVF in- 
creases with cluster velocity dispersion, while the slope is 
constant. The average slope of the observed cluster GCVF 
is somewhat flatter (/3 ~ —2.5) than that computed from 
ACDM cluster simulations (/3 — —3.4). However, the errors 
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Figure 8. GCVFs for the SDSS cluster sample and the field. The 
thin lines represent the GCVFs for SDSS clusters, normalized 
at Iog]^g(K;/ ltms~^) fs! 2.3. The solid bold line is the similarly 
norma lized observed field GCVF as determined bv lGonzalez et all 
i200(J) . The points with error bars are the field GCVF determined 
from SDSS data as described in i|4.5l Because of the arbitrary 
normalization, this plot is suitable only for assessing differences 
in the shape of the GCVF in clusters and the field. Although the 
cluster GCVFs display a range of slopes, they are steeper than 
the field for K ^ 200 and more shallow for Vc > 200. The field 
GCVFs are consistent with one another within the error bars at 
Vc < 200 km s~i. The field GCVF w e determined lies above that 
measured by iGonzalez et all i2OO0h at Vc > 200kms-l because 
we took into account the scatter around the relations used to 
compute Vc for early type galaxies, which dominate the GCVF 
at large Vc- 

on these average slopes are large, and overlap. If the differ- 
ence between the simulated and the observed clusters is real, 
it is most likely due to the fact that the simulations do not 
include the effects of baryons on galaxy rotation curves. 

4. We estimated the effects of baryons on the simulated 
GCVFs. Baryons can flatten the simulated GCVF enough to 
bring it into agreement with the GCVFs of observed clus- 
ters with similar velocity dispersions. This agreement be- 
tween the amount of cluster substructure predicted by N- 
body simulations an d that observed in r eal clusters has been 
found previously bv lMoore et al.l 1^9^ for an SCDM sim- 
ulation of a Virgo-like cluster. Our approach is different in 
that we treat early and late types distinctly, as well as ac- 
count for the effects of baryons on galaxy rotation curves. 
In addition, our simulated cluster has evolved in the concor- 
dance cosmology. 
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APPENDIX A: POPULATING SUBHALOS 
WITH GALAXIES 

In order to estimate the effects of baryons on galaxy circular 
velocity profiles, we must know the circular velocity profiles 
of both the dark matter halo and the embedded collapsed 
baryons. A full description of the baryonic mass distribu- 
tion requires knowledge not only of the density profile shape 
(galax;y type), but also a normalization (total galaxy mass) 
and a scalelength (galaxy size). In the remainder of this 
section, we outline our method for assigning galaxy types, 
masses, and sizes to the population of subhalos identified 
in the zero-redshift output of our highest resolution cluster 
simulation. Once the dark matter and baryon density pro- 
files are fully specified, it is straightforward to compute the 
final rotation curve due to both components. 

AO.l Assigning Galaxy Types 

Although galaxies exhibit a range of bulge-to-disk ratios, 
for simplicity we assign either a pure spheroidal or a pure 
disk baryonic component to each subhalo in the HR Virgo 
simulation. Sphoroidals are described by a Hernquist density 
profile (Hcrnauist 1990): 

2^ r (r + i?.ph)^' 
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where M^ph is the total mass of the spheroid, and the scale- 
length i?aph is related to the projected half-light radius -Ro of 
a deVaucoleurs surface brightness profile by i?sph = 0.551-Ro- 
Late type galaxies are assumed to be infinitely thin, and 
have surface density profiles of the form 



E(r) = 



Md 
2^ Hi' 



r/Rd 



(A2) 



where Md is the total mass of the disk and -Rd is the disk 
scale length. 



Md = 27rSoi?d- 



(A3) 



Observationally, galaxy type is a function of circu- 
lar velocity, with early type galaxies dominating at high 
Vr- JShimasakulll993l: iKochanek fc Whitell20oJ : ISheth et alJ 
I2OO3I) . Only the peak circular velocity of the dark matter 
component, Vpk, can be measured for each subhalo. We 
therefore assume that if early type galaxies dominate at high 
Vc, they also dominate at high Vp^. We randomly assign each 
subhalo the designation of early or late type based upon the 
early type fraction as a function of Vpk, /cariy(H>k)- For sim- 
plicity we adopt a step function for /early (V^k). 



/early (Vpk) 



Vpk < Vbroak 
Vpk S5 Vbroak 



(A4) 



where 125 kms"^ < Vbroak < 150 kms"'^. Alternative func- 
tional forms which provide smoother transitions in the run 
of galaxy type with V^k do not produce significantly differ- 
ent results. Subtantially different values of Vbroak produce 
dis agreement with the jo int distribution of Ro and ao found 
bv lBernardi et alJ (|20o3) (see iMli . 



AO. 2 Assigning Galaxy Masses 

In the previous subsection, we described our procedure for 
assigning either a Hernquist or exponential density profile to 
each subhalo. Each of these two profiles is parametrized by 
both a mass in collapsed baryons (A/sph or Md) and a scale 
length (i?sph or -Rd). We set the baryonic mass embedded 
within each subhalo equal to a fraction of the virial mass of 
that subhalo at the time that the baryons collapsed. 

The first step towards determining the mass in cooled 
baryons embedded within each subhalo is to determine the 
virial mass of each subhalo. The virial radius of a simulated 
dark matter halo is usually taken to be the -R200 , the radius 
interior to which the average density is equal to 200 times the 
critical density. The virial mass is then estimated as M200, 
the mass interior to -R2oo- For subhalos, -R200, and therefore 
M200 is very difficult to measure. In addition, the baryons 
likely collapsed prior to the subhalo falling into the cluster 
and suffering the dynamical effects therein. We are there- 
fore interested in the value of M200 for each subhalo before 
cluster incorporation. Although the cluster environment can 
lead to significant mass stripping, the peak circular veloci- 
ties o f halos (Vpk) remain relatively constant fHavashi et alj 
|20o3). We can therefore use the values of Vpk measured in 
the redshift zero simulation output to infer M200 using a cor- 
relatio n between the two found bv lNavarro. Frenk fc White! 



log 



M2. 
10^0 Mp 







3.23 log 



V kms-iy 



5.31. 



(A5) 



Once M200 has been determined for each subhalo, 
the second step is to specify the fraction, /cool of M200 
that winds up as cooled baryons. For late-type galax- 
ies, w e adopt the fitting function found by Gonzalez et alJ 
1I2OOOI) in an analys i s of the semi-analytic models of 
ISomerville fc Primackl (Il999l) : 

/eool^O.l(^i^), (A6) 

where x = Vpk/ (200 kms~^). The above is valid for late 
type galaxies with Vpk ~ 60 - 350 kms~^. For subhalos with 
Vpk < 60kms~\ we set /oooi = 0.02. This minimum value of 
/cool does not affect our comparison to observations, as sub- 
halos with Vpk ^ 60 km s~^ are unlikely to have Vc > Vcompi. 
For early type galaxies, we explore /cool = 0.025, 0.05, 0.1. 



A 0.3 Early Type Galaxy Sizes 

In the previous two sections, we presented simple prescrip- 
tions for assigning a density profile shape and normaliza- 
tion to each subhalo within our highest resolution cluster 
simulation. Before the effects of baryons on galaxy rotation 
curves can be explored, we must assign scale lengths to the 
baryonic component of each subhalo. One projection of the 
Fundamental Plane tells us that the line of sight velocity 
dispersions, aios, or early type galaxies are correlated with 
their half-light radii, -Re, which are in turn simply related to 
the scale length of a Hernquist density profile. Under the as- 
sumption that (Jios is also correlated with the peak circular 
velocities Vpk of the dark matter halos in which they form, 
we can relate -Rsph to V^pk. 

In order to assign values of -Rsph to each subhalo, we 
first assign values of Rq. Using the full covariance matrix 
from B03, we draw one (aios, -Ro) pair for each subhalo se- 
lected through the procedure described in ijAO.ll to host an 
early type galaxy. Each value of -Ro is assigned a rank, with 
the largest -Ro receiving a rank of 1. Ranks are similarly as- 
signed for each drawn crios and for each subhalo based opon 
its measured Vpk. Although we do not know the relationship 
between crioa and Vpk, we expect both to correlate with -Rq. 
Thus, we assign each (crios, -Ro) pair to a subhalo that has 
the same rank as that of the crios in that pair. This method 
of assigning sizes to subhalos preserves the scatter in the 
correlation between ctios and Ro, but assumes a perfect cor- 
relation between crios and Vpk. 

Once values of -Ro have been assigned to each subhalo, 
-Rsph is determined through -Rsph ~ 0.551-Ro. 



AO. 4 Late Type Galaxy Sizes 

Our statistical assignment of early type galaxy sizes to sub- 
halos was driven by observed correlations. For disk galax- 
ies, if baryons initially share the spatial and angular mo- 
mentum distribution of their host dark matter halo, and 
conserve angular momentum as they collapse, the result- 
ing scalel ength is related to the properties of the dark mat - 
ter halo JPalcanton et al.lll997l: iMo. Mao fc White! Il998f) . 
We use this idea to assign a scale lengt h, -Rd, to each sub- 
halo chosen to host a late type galaxies. !Mo. Mao fc Whit3 
(11998) show that 



-Rd = 



jd , r-l/2j. 

: Ar200/c .Tr, 



(A7) 



The Cluster Galaxy Circular Velocity Function 



where jd is the fraction of the total halo angular momentum 
contained in the disk; /cool is, as before, the fraction of the 
dark matter mass contained in cooling baryons; A is the spin 
parameter of the dark matter halo; and r2oo is the virial 
radius. Additionally, 




where c is the halo concentration and 



' f \ \ (-0-06+2. 71/„„i+0.0047/,„„i/{idA)) 

Z£2£iA] (A9) 
X (l - 0.019c + 0.00025c^ + 0.52/c) (AlO) 



We h ave specified f cooi in ilA0.2l and we assume that jd = 
f^^^i . iBrillock etlil] l^2001^ analyze a 60 Mpc simulation 
carried out with the adaptive refinement tree N-body code, 
and found that the average halo spin parameter, A — 0.04, 
independent of Vpk. They also find the following relation for 
halo concentration: 



14,/ Yp^^ (All) 

200kms-i ^ ^ 

Finally, r2oo is the radius within which the average halo 
density is 200 times the critical density. It can be calculated 
straightforwardly from M200, the mass enclosed within this 
radius: 

where pcrit.o ~ 136MQkpc~^ is the critical density of the 
universe and M200 is determined from Equation i|A5l 



